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We study the properties of a simple lattice model of repulsive particles diffus- 
ing in a pinning landscape. The behaviour of the model is very similar to the 
observed physics of vortices in superconductors. We compare and discuss the 
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, equilibrium phase diagram, creep dynamics, the Bean critical state profiles, 

O I hysteresis of magnetisation loops (including the second peak feature), and, in 

o 



particular, "aging" in relaxations. 

Important dynamical phenomena ranging from slow relaxations or hysteresis, to the 
anomalous "second peak" in magnetisation loops, are found in vortex physics of many differ- 
ent superconductors within a broad range of material parameters. This observation suggests 



^ ■ that some basic general mechanisms are responsible for the observed phenomenology 

and that schematic models from statistical mechanics can be successfully used to describe 
vortex matter 

We consider here a simple statistical mechanics model that appears to reproduce a very 
wide range of properties of vortices, ranging from dynamical behaviours to phase transitions. 
The model is an extension of a Multiple Occupancy cellular-automaton-like Model recently 
introduced by Bassler and Paczuski (BP) to study vortex dynamics at coarse grained level. 
We introduce the vortex Hamiltonian in order to be able to consider non-zero temperature 
effects in a consistent way and study them by Monte Carlo (MC) and replica theory methods. 
Our extension of the BP model also limits the occupancy of the individual lattice sites to 
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correctly take into account the finiteness of the upper critical field. This point is of crucial 
importance for the phenomenological predictions of the model. This leads us to a Restricted 
Occupancy Model (ROM). 

We find that even the two dimensional version of the model is able to qualitatively 
reproduce many features similar to those observed in real superconducting samples, including 
a reentrant equilibrium phase diagram, creep dynamics, hysteresis of magnetisation loops, 
"second peak", and others. Here, in particular, we describe its off equilibrium magnetic 
properties. The model, simple and thus tractable, nevertheless appears to capture significant 
aspects of the essential physics and help to establish a simple unified reference frame. 

The model - A detailed description of the interaction potential, U{r), between vortices 
depends on the considered region in the temperature-magnetic field (T — H) plane. For 
instance, at low field the London approximation can be used to derive two body potentials 
0], whereas at elevated fields other approximations, such as the lowest Landau level approx- 
imation, may become relevant (see eg. Like in the BP model we consider here a coarse 
grained lattice version of an interacting vortex system, with a coarse graining length scale, 
/o, of the order of the natural screening length of the problem (typically, the magnetic pen- 
etration length A). After coarse graining, the original interaction potential, U, is reduced to 
an effective Hamiltonian coupling A. In this way a drastic reduction of degrees of freedom is 
accomplished and the resulting schematic effective model can be more easily dealt with. The 
price to pay is the the loss of information on scales smaller than Iq. However, some general 
features of the system behaviour can survive at the level of the coarse grained description. In 
the above perspective, below we only consider the essential properties of vortices interaction, 
i.e., a mutual repulsion amongst vortices together with a spatially inhomogeneous pinning 
interaction. The present description can be, of course, refined by reducing the value of Iq. 
We consider the Hamiltonian: 

ij i i 

In eq.(0), rij G {0, Nc2} is an integer occupancy variable equal to the number of particles 
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on site i. The parameter Nc2 importantly bounds the particle density per site below a critical 
value and represents the upper critical field Bc2 in type II superconductors. Particles also 
have a "charge" Si = ±1 and neighbouring particles with opposite "charge" annihilate. 
The first term in eq.(0) represents the repulsion between the particles 0. Since the coarse 
graining length is taken to be of order A we choose a finite range potential: An = Aq] 
Aij = Ai if i and j are nearest neighbours; Aij = for all others couples of sites. The 
second term in eq.(|l]) just normalises the particle self-interaction energy. The third term 
corresponds to a random pinning potential, with a given distribution P{A^), acting on a 
fraction p of lattice sites (below we use p = 1/2). For simplicity we choose a delta-distributed 
random pinning: P{A^) = (1 — p)6{AP) + p6{AP — Aq). To control the overall system 
"charge density" we can add a chemical potential term — /i J2i Si to the above Hamiltonian 
{Si = SiUi). The parameters entering the model can be qualitatively related to material 
parameters of superconductors. The inter-vortex coupling Aq sets the energy scale. The 
ratio = Ai/Aq can be related to the Ginzburg-Landau parameter k = A/.^ ||^ and, in 
general, is expected to be an increasing function of /t. The last parameter A^ is a fraction 

of Aq. 

To understand the equilibrium properties of the ROM model we briefiy consider its replica 
mean field theory (MF). In this approximation the equilibrium phase diagram in the plane 
{H*,T*) (where T* = ksT/Ai and H* = ^i/keT) can be analytically deah with (see Fig.p. 
In absence of disorder it clearly shows a reentrant phase transition from a high temperature 
low density fiuid phase to an ordered phase, in analogy to predictions in superconductors 

For moderate values of the pinning energy (Aq < Ai), a second order transition still 
takes place, which at sufficiently strong pinning is expected to become a "glassy" transition, 
as is seen in Random Field Ising Models [0. For the 2D lattice we consider below (in limit 

^ 0), a numerical investigation is consistent with a first order transition. In MF, the 
extension of the low T phase shrinks by increasing A^ (i.e., the highest critical temperature, 
T^, decreases) and the higher is n* the smaller the reentrant region. These findings are in 
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agreement with experimental results on vortex phase diagrams (see Ref. [0 or, for instance, 
2H-NbSe2 superconductors from Ref. [p!T[] ). 

We now go beyond MF theory and discuss the dynamical behaviour of the model. We 
performed MC simulations on a 2D square lattice system (we use typically = 32^) de- 
scribed by eq.([l|). The system is periodic in the y-direction and has the two opposite edges 
in the ^/-direction in contact with a reservoir of particles. The reservoir is described by Ti 
with = \fi and kept at a given density Next- Particles undergo diffusive dynamics and 
are introduced and escape the system only trough the reservoir. The parameters of our 
simulations are usually Aq = 1.0; Aq = 0.3; Nc2 = 27. We have sampled several values of 

K* e [0,0.3]. 

We are interested in the dynamical properties of the system in the low T* region of the 
above phase diagram. Here the 2D model has interesting magnetic hysteretic behaviours. 
In our MC simulations we ramp Next (starting from zero and later back to zero) at a given 
rate 7 = ANq/t and record the magnetisation, M = Ni^ — Next {Ni^ = ^i^i) / is 
the "charge" density inside the system) as a function of Next- Such a ramping induces a 
Bean-like profile in our lattice (inset of Fig.|I|) with a structure similar to some experimental 
data (see, for instance, |T^). 

At low temperatures (T < 5 |]l3l), a pronounced hysteretic magnetisation loop is seen 
(see Fig. and when k* is high enough (k* > 0.25) a definite second peak appears in M. In 
the present case the origin of the second peak is very simple. At high density and k*, groups 
of vortices, frustrated in minimising their repulsive interaction energy, are forced to cluster 
together forming macroscopically extended energetic barriers which cage other diffusing 
vortices. In Fig.^we plot the average energy barrier, AE{Next), a particle meets during the 
same runs for M shown in Fig]^. A "trial" vortex approaching groups of clustered vortices 
has to pass over these barriers to move further. This dynamically generates the second peak. 
The final decrease in M at high Next is, here, due to a "softening" of these barriers caused 
by saturation effects related to the finite value of Ne2- The first peak in the magnetisation 
stems from the fact that density variations in the reservoir are only slowly transmitted in 
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the system when it is in the low density "fluid" phase. The second peak and hysteretic loops 
at moderate-high n* are also present when — > (Aq also determines the difference in the 
amplitude of |M| between the increasing and decreasing ramps). Very similar magnetisation 
data are observed in a number of different superconductors from intermediate to high k 
values (see, for instance, ref.s in P JTT|JT^ ). 

The actual shape of loops strongly depends on the parameters of the dynamics (and 
system size). In particular, the sweep rate of the external field, 7, is very important. As 
soon as the inverse sweep rate is smaller than the characteristic relaxation time (which 
can be extremely long, inaccessible on usual observation time scales, see below) strong off- 
equilibrium effects are present, such as metastability or "memory" and "aging" [^,|Tl],|l4 



As a first example of these facts, we show in the right inset of Fig.|^ the dependence of the 
second peak height, Mp, on 7. At low temperatures T < 1 and not too low /t* (/t* > 0.28), 
Mp is approximately logarithmically dependent on 7 over several decades: 

Mp(7) ~Mo + AMln(7) (2) 

Such a behaviour gets closer to a power law Mp ~ 7^ at lower k* (for instance, x ~ 1/2 
at K* = 0.26). Eventually, when 7 is smaller than a characteristic threshold, 7^, Mp expo- 
nentially saturates to its asymptotic value (this usually is orders of magnitude smaller than 
Mp(7) at high 7). Interestingly these findings are also very close to what is experimentally 
observed in superconductors an example from an YBCO sample (from |15|) is given 

in the left inset of Fig.^. The threshold, 74, is strongly dependent on the system density 
Nin (and system size) and is a rapidly decreasing function of k*; for instance at T = 0.3, 
for Next = Np {Np is the location of Mp ), 7fU.=o.26 ^ 4.5 • 10"^ but 7fU.=o.28 < 10"^ 
•y^^iT, Nin, K* , L) is a measure of the system characteristic equilibration times (which can 
be huge). 

Seemingly a dynamical phenomenon, in the ROM model the second peak is related to 
a true transition: in the 7 — limit, its location, Np, is associated with a sharp jump in 
Meq = lim^^o^(7), where its fluctuations increase with system size (see inset of Fig.^l). 
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These findings are consistent with experiments (for instance, see Ref. [|TT|) and to some 



extents reconcile opposite descriptions ("static" v.s. "dynamic") of the phenomenon. 

It is also interesting to consider the "creep rate" Q = ^j^, which is often associated 
to a measure of the intrinsic energy barriers in the creep process ||^. Experimentally, Q 
is a non trivial function of the magnetic field (see for instance We find that, due 

to the very long relaxation times, Q{N(,xt) is in itself a (slowly varying) function of 7, up 
to when 7 is smaller than the smallest 7t. In Fig.§ we show how Q depends on 7 in the 
ROM model: for T = 0.3 and K* = 0.28 we plot as a function of Next the average of Q over 
two different 7 intervals 7 G [5 ■ 10"^, 10"^] (filled circles) and 7 G [5 ■ 10~^, 5 ■ 10"^] (filled 
squares). The difference between the two is apparent. We note a remarkable correspondence 
with experimental data for YBCO, shown for the quoted sample in the inset of Fig.^. 

The presence of the above "memory" effects indicate that the system, on the observed 
time scales, can be well off equilibrium. To reveal the underlying non-stationarity of the 
dynamics we consider two times correlation functions and, at a given N^xt, we record {t > t^) 



C(t,U = ([iV,„(t)-iV,„(tJ]2) . (3) 

Fig.|| clearly shows that C(t, t^) exhibits strong "aging" : it explicitly depends on both times, 
contrary to situations close to equilibrium, where C is a function of the times difference t—ty^. 
In particular at high and low T (where relaxation times are very high) C{t, t^) can be well 
fitted by a generalisation of a known interpolation formula, often experimentally used [|l|, 
which now depends on the waiting time, t^: C{t,tw) — C^o \^ ~ -'- + fr^'^ ( t^+to ) ^^}" 
We found that to take yU ^ 1 is consistent with our data. Notice the presence of scaling 
properties: for not too short times C is a function of only the ratio t/t^\ C{t, t^^) ~ S{t/tj^). 
This is a fact in agreement with general scaling in off equilibrium dynamics (see Ref. |]I8[) 



and in strong analogy with other systems (from glass formers to granular media |T7|-pO[|). 
Experimental measurements of C{t,tw) would be very valuable. 

In conclusion, in the context of a simple tractable model we depicted a panorama of 
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magnetic properties of vortices in type-II superconductors. Even the 2D version of the 
model has many interesting features in correspondences with experimental results and allows 
clear predictions on the nature vortex dynamics. The origin of the slow off equilibrium 
relaxation (observed at low T) is caused by the presence of very high free energy barriers 
self generated by the strong repulsive interaction between particles at high densities (for k* 
above a threshold). In this respect the pinning potential plays a minor role. For instance, 
the presence of a "second peak" in M, is also observed in the hmit — > 0. The second 
peak doesn't mark the transition to a "glassy" phase, but is also present in such a case. 
sets the position and amplitude of the reentrant order-disorder transition line, which is in 
turn distinct from the second peak locations. 

At low temperatures on typical observation time scales, the system is strongly off equi- 
librium. This is most clearly seen from "aging" found in two-times correlation functions. 
These obey scaling properties of purely dynamical origin Experimental check of these results 
would be extremely important to elucidate the true nature of vortex dynamics. 
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FIGURES 
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FIG. 1. Main frame The mean field phase diagram of the ROM model in the small pinning 
strength regime {A^ < Aq), in the plane {H*,T*) (here H* = n/ksT T* = T/Ai are the dimen- 
sionless chemical potential and temperature), for k* = 10 and Aq = 0.0; 0.5; 0.75 (res. full, dotted 
and dashed lines) and k* = 3.3 and Aq = 0.0 (long dashed line). Inset The magnetisation profile, 
M{x), as a function of the transversal spatial coordinate x/L (L is the system linear size), recorded 
while ramping the external field, A^ext (for the shown values), in the 2D ROM model (k* = 0.26, 
T = 0.3, 7 = 1.1 10~^). Notice the change in shapes for Next smaller or larger than Np 2± 13.5 
(filled v.s. empty symbols). 
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FIG. 2. Main frame The magnetisation, M, as a function of the applied field density, N^xt^ 
in the 2D ROM model for k* = 0.1, 0.26 at T = 0.3 and the shown sweep rates 7. The locations of 
the peak, different in the increasing and decreasing branches, depend on 7 and approach the same 
value in the limit 7—^0. Inset The equilibrium value of M (i.e., when 7—^0) for k* = 0.26 at 
T = 0.3. 
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FIG. 3. Main frame The average energy barrier, AE, a particle meets during diffusion in 
the same ROM lattices of Fig.§. Inset right The second magnetisation peak height in the ROM 
model as a function of the sweep rate 7 for k* = Ai/Aq = 0.28 at T = 0.3. Inset left The second 
magnetisation peak, Mp (Am^ x 10'^), in a single crystal of YBa2Cu408 at temperature 20K as a 
function of the sweep rate, 7 (mT/sec) (from Ref.[15]). 
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FIG. 4. Main frame For the same ROM models of the right inset of Fig]^ with k* = 0.28, 
the "creep rate" Q = averaged over the intervals 7 G [5 • 10~^,10^^] (filled circles) and 

7 G [5 • 10^^,5 • 10^^] (filled squares), is plotted as a function of N^xt- For comparison, a cor- 
responding magnetisation loop (M — > M/4 to have clear scale on y-axis) is also shown (dotted 
line). Inset The creep rate (for 7 G [10, 100]mT/sec, full line) and magnetisation loop (4- 10^ Am^, 
dotted line) as a function of the external magnetic field /xo-ff (T) in the same YBa2Cu408 sample 
of Fig.|. 
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FIG. 5. Time relaxation of the two-times vortex-density correlation function, C{t,tyj), in the 
2D ROM model, recorded at T = 0.1 {k* = 0.28) for the shown tw at Next = 4, 16 (resp. inset, 
main frame). Continuous lines are logarithmic fits. 



